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ABSTRACT 

The continuous need towards improving the capacity of magnetic storage devices requires 
materials with strong perpendicular magnetic anisotropy. FePd, CoPd and their Co(Fe)Pt 
counterparts very attractive candidate for such purposes. The magnetic properties of these 
films are largely dependent on the orientation and local distribution of the Llo FePd phase 
fraction which are mainly controlled by diffusion processes and involve diffusion paths of a 
few angstroms. Highly ordered as well as disordered epitaxial isotope-periodic 57 FePd/ nat FePd 
(001) thin films were prepared by molecular-beam epitaxy on MgO(OOl) substrate. Short 
range diffusion of different phases in FePd thin film induced by room temperature 130 keV 
He + irradiation was investigated at fluences up to 30xl0 15 ions/cm 2 . Conversion electron 
Mossbauer spectroscopy and synchrotron Mossbauer reflectivity experiments showed that the 
inter-atomic diffusion across 57 FePd/ nat FePd interface occur mainly via the iron rich regions. 
The ratio of the diffusion length in the Llo, fee and iron rich structure are 1:1.4:5.6 
respectively. Assuming that the diffusion coefficient in the fee phase is between the diffusion 
coefficient of the Llo phase in the crystallographic c direction and perpendicular to it, the 
diffusion coefficient in the c-direction of the Ll phase is found more than 1.9 times lower 
than in the a-b plane. 

PACS 75.70.-i, 75.30.Gw, 75.50.Ss, 61.80.-x, 74.25.Ha 



I INTRODUCTION 

Due to their high perpendicular magnetic anisotropy, (PMA), Llo (CuAu(I)-type) ordered 
FePd as well as FePt and CoPt are candidate materials for future ultra-high density magnetic 
recording [1,2,3]. PMA has been reported to originate in these materials from growth-induced 
heterogeneity, compressive strain, short-range order driven segregation or magnetically 
induced phase separation, but was not found to be linked to the appearance of any long range 
ordered phase. At room temperature, Fe!_ x Pd x exhibits equilibrium Ll structure in the 
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composition range 0.5 < x < 0.6, with alternating Fe and Pd planes along the [001] direction. 
In the Llo FePd structure the short and long lattice parameter ratio (former along the [001] 
axis) varies between 0.96 and 0.97. The anisotropy energy is in range of 10 erg/cm [4,5,6] 
which is of primary interest for magnetic storage applications since it allows nano-size 
magnetic volumes to remain ordered against thermal agitation. The ordered Llo to disordered 
fee phase transition is driven by short range diffusion. 

Atomic diffusion governs the structural relaxation and associated changes in the physical and 
magnetic properties in these alloys. Depth profiling techniques like radioactive tracer and 
secondary ion mass spectrometry are most widely used for such studies. However, the depth 
resolution of these techniques is often limited to a few houndred angstroms. This is a severe 
limitation, since in metastable systems such large diffusion lengths may be difficult to achieve 
without fundamental structural transformations. Dynamic techniques like Mossbauer 
spectroscopy or nuclear forward scattering of synchrotron radiation, on the other hand, are 
limited by the data acquisition time and the time window over which they are sensitive to 
atomic diffusion. For example, diffusivities accessible using 57 Fe Mossbauer spectroscopy are 
in the range of 10" to 10" m /s [7] which is not enough to follow diffusion in atomic range. 
In studies of chemically modulated multilayers x-ray diffraction (XRD) and x-ray reflection 
(XRR) methods are commonly used. By XRR, it is possible to study diffusion lengths below 
the detection limit achievable by the sectioning and profiling techniques [8,9,10,11,12]. 
Despite the difficulties that arise due to the concentration dependence of the diffusion 
coefficient, these techniques are widely used to study interdiffusion in compositionally 
heterogeneous multilayers, it is not possible to follow self-diffusion by these techniques in 
chemically homogenous materials since the solely electronic interaction with the x-ray 
photons produces no contrast between the adjacent layers. 

To investigate self diffusion, an isotope sensitive technique is required. Conventional x-ray 
techniques don't have this feature but neutrons may have significant difference in scattering 
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cross section for the different isotopes of the same element. Neutron reflectivity is a suitable 
non-destructive method to study self-diffusion in a chemically homogenous isotope multilayer 
with diffusion lengths of the order of a few angstr0ms [13,14]. A recently emerged similar 
isotope sensitive non destructive method for studying self diffusion is Synchrotron Mossbauer 
Reflectivity (SMR) [15]. Due to the nuclear (hyperfine) interaction between the atomic nuclei 
and the highly monochromatized synchrotron radiation (tuned to the transition energy of the 
resonant isotope, hereafter Fe), Mossbauer isotopes have very large scattering amplitude as 
compared to that of non-Mossbauer-isotopes of the same material providing an isotope 
contrast between the adjacent layers of the sample. This results in Bragg peaks in the 
reflectivity pattern, related to the period of the periodic multilayer sample. The sensitivity of 
these methods are suitable for the layer thickness used in our case. From the variation of the 
Bragg peak intensity, valuable information can be gathered on the change of the isotope depth 
profile, and consequently on the interlayer diffusion. However, isotope sensitive neutron and 
Mossbauer reflectometry has been already used by different groups [16,17] they assumed that 
the sample structure remained unchanged during the heat treatment thus neglecting the effect 
of phase evolution. 

Low energy(~130keV) He + implantation features low collision cross section with the 
intention to avoid defect interactions and small energy transfer to minimize recoil 
displacement. However, the beam energy is sufficient for the ions to pass through the layer 
and stop deep (-1.5 urn) in the substrate, leaving a rather homogeneous defect distribution in 
the (65 nm thick) film to relax. According to SRIM [18] simulations each incoming ion 
generates ~6 displacements on average, which mainly generate vacancies and equal number 
of interstitials According to these simulations, the probability of sputtering by these He + ions 
is negligible (< 10" 5 ). This kind of low energy ion implantation doesn't induce collision 
cascades in the layer and results similar effect as heat treatment. 
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In this paper we use SMR method to study the self-diffusion of Fe in an isotope periodic 
natFePd/ 57 FePd multilayer film following room temperature low energy He + irradiation of 
various doses. Isotope layers with originally sharp interfaces were used to avoid chemical 
effects of diffusion. X-ray diffraction (XRD), conversion electron Mossbauer spectroscopy 
(CEMS), high resolution transmission electron microscopy (HRTEM), completed by selected 
area electron diffraction (SAED) and SMR were used to characterize the samples and to 
follow the irradiation induced processes in the samples. 

II THEORETICAL CONSIDERATION 

SMR experiments can give information on the atomic movement below the nanometer scale, 
which is indispensable to understand self-diffusion-controlled processes. The Fe 
concentration in one isotope layer of nat Fe47Pd53/ 57 Fe47Pd53 system can be expressed in a 
Fourier series [11,19] 

C(z) = lC,«ap(ft,z) (1) 

n 

Where k n =2mi/A with A being the periodicity of the multilayer and z the depth of the sample. 
When atomic diffusion takes place in the sample by annealing or irradiation the magnitude of 
the Fourier component C„ decays with time. If the atomic migration by annealing results 
diffusion in the sample then we can use the ID Fick equation. Substituting (1) into this 
equation one obtain the following equation. 

C n =C 0n exp[-k 2 n Dt] (2) 

Where D and t are the diffusivity at a given temperature, and the retention time, respectively. 
If the atomic migration is diffusion then D depends from temperature Arrhenius like. If the 
atomic migration is irradiation induced then D depends from the energy loss of the ion beam 
for a given depth. 
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T. Mizoguchi et al.[ ] derived a formula, which relates the intensity of the n-th order Bragg 
peak in the reflectivity pattern at t = (Jo) to that at time t, according to the following 
expression: 



log 



V 'o J 



^DCDt, (3) 

A 



where log stands for the natural logarithm. The corresponding diffusion length due to 
thermally induced diffusion is: 

L d = 4lnD(T)t (4) 

In the following we show that low energy ion-irradiation can be treated in a very similar way. 
Materials exposed to ion irradiation exhibit significant atomic rearrangement, so-called ion- 
beam mixing. Several processes are responsible for the ion mixing effect, all of which are 
initiated by the interaction of energetic ions with the solid. The relative significance of the 
ballistic, recoil and cascade effects can be altered by changing the mass and/or energy of the 
ions impinging on the sample. Increasing the mass of the ions increases the amount of energy 
deposited in nuclear collisions per unit length traveled by the ion. Consequently, the amount 
of mixing, Q, at the interface of two layers can be expressed as [21]: 



Qx 







1/2 
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Where (dE/dx) n is the nuclear stopping power and <f> represents the ion dose that have passed 
through the interface. Electronic interaction between the beam and the atom results only 
ionization and causes no mixing hence it doesn't play role in our case. Since the dose rate, in 
ion/cm /s, can be considered constant during an ion mixing experiment, the ion dose is 

i yo 

proportional to time, leading to the observed behavior that mixing is proportional to (dose) 
[22]. This implies that ion mixing is also proportional to the square root of ion mixing time. 
One can compare this proportionality to that observed for a reaction layer formed between 
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two materials by thermally activated interatomic diffusion [23]. The width of the reacted layer 
W, found to behave as follows: 

W oc 4lEt = L d (6) 

where D is of the form of an interdiffusion coefficient. Because of these similarities many 
ion mixing models are based on diffusion and interdiffusion analogues with mixing described 
in terms of an effective mixing diffusion coefficient, D. 

Ill EXPERIMENT 

Highly ordered (majority Llo) and also disordered (fee) isotope-periodic Fe47Pds3 were 
grown on 20x20x2 mm 3 MgO(OOl) substrates by the method of molecular-beam epitaxy 
(MBE). In the case of ordered sample the substrate was held at 350°C and in the case of 
disordered sample the substrate was at room temperature during growth. In order to obtain 
epitaxial growth of Pd, a seed layer of 3 nm Cr was evaporated onto the MgO(OOl) surface at 

o 

a rate of 0.07 A/s, followed by the growth of a Pd buffer layer of 15 nm thickness at a rate of 
0.2 A/s. Then the bi-layer sequence of 3 nm nat Fe47Pd 53 and 2 nm 57 Fe 4 7Pd 53 was repeated ten 
times. All nat Fe47Pd 53 and 57 Fe 4 7Pd 53 layers were prepared by co-evaporation of Fe and Pd at a 
rate of 0.0485 A/s for both 57 Fe and nat Fe and 0.0682 A/s for Pd. To avoid oxidation of Fe a 1 
nm Pd layer was grown on top of the sample. The 57 Fe was evaporated from Knudsen-cell (at 
1410 °C), all other layers were deposited using electron-beam evaporation. The base pressure 
was 2xl0" 10 mbar in the MBE chamber, which raised to 2.8xl0" 9 mbar during the growth. It is 
very important to achieve identical layer composition in the nat FePd and 57 FePd layers, 
therefore the so-called MBE tooling factors were carefully calibrated using Fe and Pd as well 
as Fe x Pd y test films, the layer thickness and the composition of which were determined by x- 
ray reflectometry and Rutherford Back-Scattering. The evaporation was controlled by two 
independent quartz thickness monitors. 



Preliminary release Not for distribution 

The samples were cut into eight equal pieces of 10x5x2 mm 3 . Before He + -ion irradiation, 
SRIM [18] simulation were performed and the energy of the ion was chosen so that the He + 
passes through the layer and stops deep in the substrate without causing collision cascades. 
The ordered samples were then irradiated by l.OxlO 14 to 1.49xl0 16 at/cm 2 and the disordered 
ones were irradiated by 3.7xl0 14 to 3xl0 16 at/cm 2 He + of energy 130 keV and one sample of 
each set were left unirradiated. 

The CEMS experiments were performed by using a Co(Rh) single line Mossbauer source 
with a home made gas-flow single-wire proportional counter operating with He gas with 4.7% 
CH4 extinction gas and at a bias voltage of 830+10 V. 

SMR measurements were performed at the ID 18 and ID22 beam line of the European 
Synchrotron Radiation Facility (ESRF, Grenoble, France) in 16 buntch mode. For sequential 

o 

monochromatization of the beam to the 14.4 keV Mossbauer transition ((/l=0. 86025 A) of 
Fe, a Si(lll), then a Si(422)/Si(12.2.2) double channel cut monochromator were used. 
Prompt (non-resonant) and 57 Fe delayed (time integrated) nuclear resonant reflectograms were 
recorded on each sample. In this method the interacting photon scatter on both the atomic 
electron and nucleus. The characteristic time of the electronic scattering is typically 10" 15 s 
from which the signal of nuclear scattering is well distinguishable since it is seven orders 

Q 

slower (-10" s). This means that the delayed signal is fully originates from the nuclear 
scattering and basically free from noise. Since the monochromatization of the synchrotron 
beam is currently a few meV and this is orders of magnitude higher than the transition energy 
of the nucleus (-100 neV) the different transitions of the nucleus are excited together and this 
results interference in the time spectrum. This time spectrum can be considered as the Fourier 
transformation of the conventional Mossbauer spectrum. The Synchrotron Mossbauer 
Reflectometry is the integration of the time spectrum at each angle which gives a very useful 
isotope sensitive method. 
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X-ray diffraction was carried out using a D8 Discover type diffractometer (equipped with a 
Gobel-mirror on the primary side) in Bragg-Brentano geometry using Cu K a radiation 
</l=1.5415A). 

IV RESULT AND DISCUSSION 

The He + irradiation not only causes interface mixing in the FePd samples, but also induces, as 
we shall see, structural transformation. Detailed CEMS, XRD and TEM evaluation of the 
ordered irradiated samples has been published elsewhere [24,25]. Our earlier CEMS 
experiment show, three distinct micro structural units, the low hyperfine (hf) field, a large hf 
field, and an intermediate hf field local environment could be identified in the ordered 
samples, which were attributed to the ordered Llo, one iron-rich phase and the disordered fee 
structural units, respectively (Figure 1 a.). The iron-rich environment consists of iron 
nanoclusters within the Llo matrix, being magnetically coupled with it. In the ordered and 
iron-rich components the hf field points out of the sample plane, while the disordered phase 
has random magnetic orientation. By increasing the fluence of the He + irradiation from zero to 
14.9x10 ions/cm , the CEMS spectral fraction of the ordered Llo phase decreases from 81% 
to 44%, while the disordered phase increases from 15% to 36%. The fraction of the partially 
ordered component also increases from 3.4% to 20%, but it reaches this value already at a 
fluence of 7.4x10 ions/cm (Figure 2 a.) 
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Figure 1 The change of the CEMS spectra and hyperfine field distribution taken after growth, medium and high 
irradiation. A) ordered sample b) disordered sample. 

The CEMS characterization of the disordered samples showed the lack of the Llo and the high 
hf field component (Figure lb.) and only the disordered fee component is present in these 
kind of samples. The irradiation does not effent this structure (Figure 2 b.). 
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Figure 2 The change of fractions of the different phases in ( a) ordered b) disordered) FePd after various 
fluences of He + irradiation as obtained from CEMS spectral intensities of the different species, a) The ordered 
Ll phase decreases from 81% to 44%, while the disordered phase increases from 15% to 36%. The fraction of 
the iron ri component increases from 3.4% to 20%. b) The fraction of the disordered fee environment increased 
from 52% to 62% and the iron region decreased from 48% to 38%. 
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In order to gather information about the distribution of the local environment obtained from 
CEMS, additional selected area electron diffraction (SAED) measurements on the ordered 
sample were performed. This method is highly suitable to distinguish separate phases since it 
works in reciprocal space and so the higher index reflections of similar structures can be well 
separated On Figure 3 SAED measurement taken on the full sample (fcc-Pd + Llo) and nano 
beam diffractin taken only on Llo-FePd are compared. 
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Figure 3 Selected area electron diffraction taken from MgO, Cr, Pd and FePd layers (a) and the corresponding 
linescans taken on Ll (204) and Pd(204) (up) and on Ll (202) and Pd(202) (down). Nano-beam diffraction 
taken from only FePd layers (b) (from an area of ~10nm in diameter) and the linescan of Ll (204) (up) and 
Ll (202) down. 

When distinct fcc-FePd phases separated by phase boundaries were present in the Ll -FePd 
structure then it could be seen on the distribution of the 204 index reflection, of the SAED 
snapshot taken from the FePd layer, Since no such distribution is present we concluded that 

o 

the size of any kind of precipitate cannot exceed the size of 5 A and the structure of the 
sample can be treated as a random alloy of separate micro structural environments. 



The evaluation of the reflectivity spectra in similar previous studies [ ! ' 17 ' 26 ] was restricted to 
follow the decay of the isotopic multilayer (ML) Bragg peak height in the time integral SMR 
reflectivity curve normalized to the intensity of the total reflection peak. However, since the 
Bragg peak shape depends on the hyperfine fields (and their distributions) of the different 
species in the layers [see also 27], it is not justified to take only the normalized Bragg peak 
amplitude into account. Moreover, the normalization - in general - is only possible to the 
intensity of the total reflection peak, as an "internal standard", but its intensity is extremely 
sensitive not only to small differences in the successive beam alignments prior to recording 
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the SMR spectra of the different samples, but also, due to the limited penetration of the 
radiation, to the absorption of the near-surface layers, which may be differently effected by 
the progress of the diffusion. Consequently, normalization of the isotope Bragg peak height to 
the total reflection peak height may become rather uncertain. Therefore the in a proper 
evaluation the entire reflectivity curve has to be fitted to a layer model, which takes the 
diffusion prehistory of the multilayer properly into account. This latter procedure was 
followed here. The evaluation of the delayed SMR spectra was performed by FitSuite [28] 
program (Figure 4). The implemented routine is organized as follows. 
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Figure 4 The change of synchrotron Mossbauer refletivity spectra and the corresponding diffusion profile a) 
ordered FePd b) disordered FePd upon He + irradiation. 

The squared diffusion lengths, Wi can be summed up for successive heat treatments or 
consecutive irradiation steps. According to TEM studies the irradiated layers do not show 
distinct boundary-separated phases, therefore the Mossbauer species represent various 
environments in an otherwise single phase sample. Therefore, in the following we use the 
approximation of a multi-component random alloy. Figure 5 a and b shows the Dt evaluated 
from fitting SMR spectra for the FePd multilayers with different initial structure. On the left 
side the originally ordered sample while on the right side the originally disordered sample is 



present. 
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Figure 5 The change of the fitted Dt for the whole system after He + irradiation. A) ordered b) disordered sample. 
The dashed line indicates the Dt for a structurally homogenous system. 

At the as deposited samples there was only negligible mixing between the adjacent isotope 
layers. By the fulence of 15x10 He /cm the Dt increased upto 9 A in the case of the 
ordered sample and reached 21 A for the disordered sample. This significant difference 
indicates that the inner structure of FePd plays major role in the scale diffusion process. One 
can also see that the two curves show different tendencies with the fluence of irradiation. It 
has been shown already that the effective mixing diffusion coefficient of a structurally 
homogenous material depends linearly on the irradiation dose [29]. By the ordered sample the 
growing tendency of the Dt exhibits a slight increase which is not the case of disordered FePd 
where it shows a constant tendency with the dose. The reason for the increasing Dt is that the 
ratio of the component which blocks the diffusion are decrasing while the ones in which the 
diffusion is faster increase after irradiation. 

If we consider the model where the D piece wise constant, thus D takes an arbitrary positive value 

D j in volumes of arbitrary size, the squared diffusion length of Fe in the multilayer, 

W oc (Dt) 112 is, like in a random alloy, a weighed sums of the individual squared diffusion 

lengths of the different Fe-environments, W, = K^Dj) 1 ' 2 , i.e. the total Dt can be written as 

follows. 



Dt = x Ll Dt LXo + x fcc Dt fcc + x P0 Dt P0 



(6) 
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Where x Ll , x fcc , x P0 and Dt Ll , Dt fcc , Dt P0 are the ratios and the effective diffusion 

parameters of the ordered Llo, disordered fee and partially ordered phases, respectively. The 
ratios as well as the effective diffusion parameters depend on the dose of irradiation. The 
fluence dependence of the ratios were determined by fitting (see Figure 2). 
With these conditions (6) has the following form: 

l </> | <t> i <i> 

Dt(<f>) = (Dt ^ +a<j>)-]x Uo (<j>)d0 + (Dt^ +b<f>)-lx fcc (<f>)d<f> + (Dt° po +c^)-^x P0 ^)d</> (7) 

r o r o r o 

By simultaneously fitting the total Dt of the ordered and disordered system obtained from 
evaluation of SMR spectra (Figure 4) with (7) we get the contribution of each particular 
environment to the diffusion of the entire multilayer (Figure 6). 
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Figure 6 Contribution of the individual Fe-environments in the FePd to the total diffusion length 

The ratio of the effective diffusion parameters in the Llo, fee and partially ordered structures 
are 1:42:64, respectively. The effective diffusion lengths (perpendicular to the sample plane) 
in the different micro structural species are significantly different. Indeed, in the ordered Ll 
phase the diffusion of iron is almost blocked (the direction perpendicular to the sample plane 
coincides with the crystallographic c axis in these epitaxial samples). This means, that atomic 
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migration across 57 FePd/ nat FePd "interface" occurs mainly via the iron rich regions. 
According to (4), the diffusion lengths are related as 1:6.4:8 in the Llo, the fee and the iron 
rich regions, respectively. Taking a plausible assumption, namely that the value of the 
diffusion coefficient in the fee phase is between the diffusion coefficient of the Llo phase in 
the crystallographic odirection and the diffusion coefficient of the Llo phase in the a-b plane, 
the diffusion coefficient in the odirection of the Llo phase is more than 1.9 times lower than 
in the a-b plane. Indeed, vacancy diffusion in the odirection requires exchange of a vacancy 
between the Fe layer and the Pd layer, which have rather different vacancy creation energies. 



V. SUMMARY 

57 

In summary, the short range diffusion of Fe in highly ordered and disordered isotope- 
periodic epitaxial Fe47Pds3 thin film after varying fluence of to 14.9x10 ion/cm in the 
ordered and 0-30 ion/cm 2 of 130 keV He + irradiation in the disordered sample was 
investigated. The local structural changes were followed by Mossbauer spectroscopy. With 
increasing fluence, in the case of the ordered FePd, the fraction of the disordered and of the 
partially ordered components increased at the expense of the ordered Llo phase. In the case of 
disordered sample the Ll structure was not present and the fee region increased at the 
expense of the iron rich region. The evaluation of the full delayed SMR intensity curves 
showed considerable interdiffusion of the isotope layers. The total diffusion length of the 
multiplayer is separated to individual diffusion lengths of the distinct microstructural species 
of the otherwise homogeneous phase. By fitting the variation of effective diffusion parameters 
using the corresponding Mossbauer spectral intensities of the individual Fe-environments as 
weights, we obtained an effective diffusion length for each Fe-environment. We find the ratio 
of the effective diffusion lengths in the Ll , fee and partially ordered structure to be 1:6.4:8, 
respectively By assuming the diffusion coefficient of the fee phase to be between that of the 



Preliminary release Not for distribution 

Llo phase in the crystallographic odirection and perpendicular to it, the diffusion coefficient 
in the odirection of the Llo phase is more than 6.4 times lower than in the x-y plane. 
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